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We study the expansion of ultracold neutral plasmas in the regime in which inelastic collisions 
are negligible. The plasma expands due to the thermal pressure of the electrons, and for an initial 
spherically symmetric Gaussian density profile, the expansion is self-similar. Measurements of the 
plasma size and ion kinetic energy using fluorescence imaging and spectroscopy show that the 
expansion follows an analytic solution of the Vlasov equations for an adiabatically expanding plasma. 



Exactly solvable problems are rare in physics and serve 
as ideal models that provide a starting point for under- 
standing more complex systems. Here, we report the ex- 
perimental realization of a laser-produced plasma whose 
dynamics can be described by an exact analytic solution 
to the Vlasov equations which are central equations 

in the kinetic theory of plasmas. Expansion into a sur- 
rounding vacuum is fundamentally important and typ- 
ically dominates the dynamics of plasmas created with 
pulsed lasers ^ , such as in experiments pursuing incrtial 
confinement fusion ^] , x-ray lasers Q , or the production 
of energetic (>MeV) ions through irradiation of solids 
d, , thin foils [i i, [m, [12, lli, Q m , clusters ^ . 

We study plasma expansion with ultracold neutral 
plasmas (UNPs) [13, [I3|, which are created by pho- 
toioinizing laser-cooled strontium atoms [l3| just above 
the ionization threshold. These systems occupy a 
new regime of neutral plasma physics, and their well- 
controlled initial conditions and relatively slow dynamics 
have enabled precise studies of strongly coupled plasma 
properties [13, 21 1 . Here we use fluorescence imaging and 
spectroscopy for the first time on these systems. These 
powerful diagnostics allows us to apply the advantages of 
UNPs to a new class of important problems. 

The investi gati on of plasma expansion dates back 
many decades 22 , 2^ . Recently, exact solutions for spa- 
tially finite plasmas expanding into vacuum were identi- 
fied for a 1-dimensional plasmas [ij and later extended to 
3-dimensions [3, [31 . This work was motivated by plas- 
mas produced with short-pulse lasers. 

The Vlasov equations, along with Poisson's equation, 
describe the evolution of the electron (a = e) and ion 
(a = i) distribution functions, /Q,(r,v). The Vlasov 
equations neglect radiative processes and coUisional phe- 
nomena such as electron-ion thermalization and three- 
body recombination (25l |. but this formalism describes 
many types plasmas and is part of the foundation of ki- 
netic theory. 

Among broad classes of general analytic solutions to 
the Vlasov equations 2j|, UNPs realize a particular so- 
lution, that even describes coUisional systems and is valid 
for a quasi-neutral plasma with spherical Gaussian dis- 
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Quasi- neutrality is defined by ~ rii, where the elec- 
tron and ion densities are na{r) = J dvfa{r,v). Ta are 
the electron and ion temperatures, and the local average 
velocity varies in space according to u(r, t) = 7(t)r. The 



temperatures must scale as a^Ta — constant [l|, which is 
expected for adiabatic cooling in a spherically symmetric 
UNP [2i,[2i]. 

Under these conditions, the plasma dynamics is given 
by r„(<) = T„(0)/(1 + tV^-elp), where the char- 
acteristic expansion time r^xp is given by t, 
v/mia(0)2/A:B[7;(0) + Ti(0)]. Also, 
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describe the evolution of the characteristic plasma size 
and ion velocity, which are important experimental ob- 
servables. We define the rms 1-dimensional ion velocity 
\/ ((v • y)^) = Vi^rmsi where y is the laser propagation 
direction, v is the total ion velocity including random 
thermal motion and expansion, and the angled brackets 
refer to an average over the distribution function. 

An intuitive explanation for the self-similar nature of 
the expansion is that thermal pressure produces an aver- 
age radial acceleration [3, [3] 
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The simplification implied by the last equality is only 
valid for a spherical Gaussian plasma, and the linearity 
in r preserves the shape of the distribution functions. 

Plasmas produced with solid targets, foils, and clus- 
ters are often quasi-neutral and well-described by the 
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Vlasov equations, and electrons cool adiabatically dur- 
ing much of their evolution, but experimental condi- 
tions studied are typically very complicated and evolve 
extremely rapidly, which frustrates detailed comparison 
between experiment and theory. Final ion kinetic en- 
ergy distributions have been shown to agree with simple 
models [l^ 16 1, but in general, these systems lack the 



respectively. 



Gaussian distribution functions necessary to realize the 
analytically describable self-similar expansion. 



For appropriate initial conditions, UNPs [17| fulfill 
the requirements for the analytic solution. UNPs have 
ion temperatures of about IK determined by disorder- 
induced heating [2^, 27 1. Electron temperatures ranging 
from 1 to 1000 K are set by the detuning, E^, of the 
ionization laser above threshold. The peak densities are 
on the order of 10^^ m~^, and the profile follows that 
of the laser-cooled atom cloud, which we adjust to have 
a spherically symmetric Gaussian. The photoionization 
pulse length (~ 10 ns) is much less than the expansion 
time scale (~ 10/zs). 

The electron distribution equilibrates locally within 
100 ns and globally within I fis after photoionization 
(iil [2^. This ensures a Gaussian electron distribution 
function at the start of the expansion. Despite these 
very rapid electron-electron collisions the corresponding 
collision integral vanishes for the spherically symmetric, 
Gaussian velocity distribution eq.(l). Hence the highly 
collisional UNP considered here provides an ideal model 
system for truly coUisionless plasmas behavior. 

Ions reach local thermal equilibrium within a few 
100 ns [2^. They do not equilibrate globally on the time 
scale of the expansion 2]| , but the ions are so cold com- 
pared to the electrons, that the lack of a global ion tem- 
perature does not cause any significant deviation from 
the exact solution. The low ion temperature also implies 
that the ions form a strongly coupled fluid [l^,!!^, which, 
however, negligibly affects the plasma expansion [2^. 

The analytic expansion solution has been discussed 
previously for UNPs 2^ 28, 13], and it has been checked 
against average terminal ion expansion velocities [3l| and 
measured electron temperatures (s^ that qualitatively 
affirm the importance of adiabatic cooling for appropri- 
ate initial conditions. The lack of spatial and temporal 
resolution, available here, however prevented conclusive 
tests of the analytic predictions. Cummings ei a/. 33L 3^ 
adapted the formalism of [2^ 28, Is^l and used light- 



scattering from a small region of the plasma to study 
the expansion of ultracold plasmas with an elongated as- 
pect ratio, but they found significant deviations from the 
predictions of the model that perhaps arose because the 
condition of spherical symmetry was not fulfilled. 

To demonstrate that ultracold neutral plasmas can re- 
alize the analytic expansion solution, we will first de- 
scribe our diagnostic and show that the plasma remains 
Gaussian during its expansion. Then we will show the 
size variation and ion velocity are given by Eqs. [5] and [3] 
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FIG. 1: Recording fluorescence of UNPs. The correlation be- 
tween position and expansion velocity (red arrows) produces 
a striped image when the Doppler-shift due to expansion ex- 
ceeds the Doppler shift associated with thermal ion velocity. 

Figure [T] shows a schematic of the fluorescence imag- 
ing experiment. A laser beam that is near resonance with 
the — transition in Sr"*" at A = 422 nm prop- 

agates along y and illuminates the plasma. Fluorescence 
in a perpendicular direction (z) is imaged with a 1:1 re- 
lay telescope onto an image-intensified CCD camera. The 
422 nm light is typically applied in a 1 fis pulse to pro- 
vide temporal resolution, and the intensity is only a few 
mW/cm^, which is low enough to avoid optical pumping 
to the metastable state. 

A general expression for the fluorescence is 
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where u is the laser frequency and ^eff — To ~t~ ^laser is 
the sum of the natural linewidth of the transition (27r x 
20 MHz) and the imaging laser hnewidth {2n x 8 MHz). 
Ti.therryii^) is thc local temperature of the ions describing 
random thermal motion, which gives rise to the Doppler 
width (T£). Due to the directed expansion velocity, the 
average resonance frequency of the transition for atoms 
at r is Doppler-shifted from the value for an ion at rest, 
vq, by i^|'a;p(r) = u(r) • y/\. The spatial variation in 

Ti^therm, which varics as [2l[, is small compared to 
the directed expansion energy, so T^z^i^erm can be taken 
as constant. 

Images can be analyzed in several different ways, which 
each provide access to different plasma properties. Sum- 
ming a series of images taken at equally spaced frequen- 
cies covering the entire ion resonance is equivalent to in- 
tegrating F{iy,x,y) over frequency. This yields a signal 



3 




2.5- 



FIG. 2: Evolution of the areal density for a plasma with 
1.6 ± 0.1 X 10* ions, 2E^/2,k^ = 48 ± 3K, and cr(0) = 
1.3 ± 0.1mm. The arial density (arbitrary units) is found 
by summing together 50 images taken at equally spaced fre- 
quencies that fully cover the ion resonance. The Gaussian 
fits (solid line) to linear cuts show that the expansion is self- 
similar. The time indicated is the evolution time since plasma 
creation. The right-hand axes show that the differences be- 
tween data and fit are small. 



proportional to the areal plasma density, J n{x,y, z)dz, 
which for a Gaussian density distribution should take the 
form nareai{x,y) = \/27rCTnoexp [-(x^ +y'^)/2a'^]. 

Figure [2] shows that Gaussian fits of the areal density 
are excellent during more than a factor of two change in 
a from the earliest times until the signal expands beyond 
the range of the imaging system. This provides direct 
confirmation of the self-similar nature of the expansion. 

There is no sign of any deviation from the gaussian 
shape at large radius. This might seem surprising be- 
cause self-similarity follows from Eq. 31 and this equa- 
tion must break down at large r where it implies unphys- 
ical accelerations arising from unphysically large electric 
fields. One would expect a breakdown where the plasma 
is not quasineutral, which should occur when the local 
Debye screening length exceeds the length scale for ion 
density variation, or a <C ^nir) = [£ofcB7e/ne(r)e^]^/^ 
for a Gaussian plasma. At the radius where a ~ Xoir), 
the accelerating electric field reaches a maximum of 
^fceTe/eAi), identical to peak accelerating fields in plas- 
mas generated from foil targets [s^. A peak in the 
electric field can lead to wave-breaking, shock waves, 
and the formation of an ion front as seen in theories 
for UNPs [2^ 0, and traditional short-pulse laser- 
produced plasmas [15, [s^, 37 1. The lack of such features 
in data such as Fig.[2]confirms the extent of quasineutral- 
ity one would expect for such a low electron temperature 
isl 30t and may also suggest that finite ion temperature 
and strong coupling effects damp wave-breaking at the 
plasma edge [26]. 

The evolution of a can be extracted from fits such as 
Fig. [21 as shown in Fig. [31 for several different plasma 
conditions. Again, the data evolve as predicted by the- 
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FIG. 3; Evolution of the plasma size, extracted from fits 
such as in Fig. [1 re(0) and cr(0) are fit, while ri(0) is 
taken from the theoretical expression for disorder-induced 
heating [20I. [27|. Uncertainties in 2Es/3kB reflect 1-standard- 
deviation calibration uncertainty in the wavelength of the 
photoionizing laser. Quoted uncertainties in Te(0) are sta- 
tistical, but there is an additional systematic uncertainty of 
a few percent arising from calibration of the imaging-system 
magnification and overlap of the plasma and fluorescence ex- 
citation laser. Statistical uncertainty in the measurement of 
a is less than the size of the plotting symbols. Initial peak 
densities for these samples are ~ 10^^ m~'^. 



ory. The expansion is sensitive to the electron temper- 
ature and initial size, and higher electron temperature 
and smaller size lead to a faster expansion because this 
increases the thermal pressure (Eq. [J]). Fits of the data 
using Eq. [2l yield values of Te (0) that agree reasonably 
well with the expected values of 2Ef./'ik-Q. This con- 
firms that on the time scale of the expansion there are no 
significant collisional or radiative processes changing the 
electron temperature such as three-body recombination 
or electron-ion thermalization, as assumed in the Vlasov 
equations. 

In order to completely characterize the plasma expan- 
sion we also measure the light-scattering resonance spec- 
trum, formed from the integrated fluorescence in each of 
a series of images taken at different frequencies, as shown 
in Fig. [31\_. Eq. [51 combined with the expansion velocity 
u(r) predicted by the Vlasov equations, implies that the 
resulting signal dx dy F{v,x,y)) should take the form 
of a Voigt profile. The rms width of the Gaussian com- 
ponent of this profile arising from Doppler broadening 
reflects both thermal ion motion and directed expansion 
and is given by ctd = Vi^rms/)^- 

Equation [31 provides an excellent fit to the data, and 
the extracted values oiTe{Q) are consistent with 2Ee/3kB 
within experimental uncertainty. The small initial offset 
of Vi,rms in each plot is due to disorder-induced heating 
of the ions within the first microsecond [23|, which lo- 
cally produces a thermal ion velocity distribution. But 
yi,rms quickly increases well above this value as the elec- 
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FIG. 4: (A) Sample spectrum and (B) evolution of the rms 
ion velocity (vi,rms). (A) Spectrum for 2Ee/SkB = 105 ± 3 K, 
peak ion density of 3 x 10^^ m~^, and plasma evolution time of 
3 ps. Fits to a Voigt profile provide Vi^rms- (B) The velocity 
evolution is fit well by Eq. |3]with Te(0) as a fit parameter, 
ct(0) determined from the cloud size measurements (Fig. [Sj, 
and Ti(0) taken from the theoretical expression for disorder- 
induced heating [13, [13]. 



tron pressure drives the plasma expansion and electron 
thermal energy is converted into directed radial ion veloc- 
ity. As fT increases and the electrons cool adiabatically, 
the acceleration decreases (Eq. [J]) and the ions eventu- 
ally reach terminal velocity when essentially all electron 
kinetic energy is transferred to the ions. 

Our measured density profiles (Fig. [3|) confirm the va- 
lidity of Eq. [3l which shows that the ion acceleration is 
sensitive at all times to the instantaneous electron tem- 
perature and width of the plasma. Agreement between 
experiment and theory for both the size evolution and the 
ion velocity dynamics hence demonstrates that Tf.{t) also 
follows the dynamics predicted by the analytical solution 
to the Vlasov equations. 

We have demonstrated a plasma in which the expan- 
sion matches an analytic solution to the Vlasov equa- 
tions [2| proposed as a basic model for the dynamics 
of quasineutral laser-produced plasmas. To realize this 
situation experimentally, it is necessary to create UNPs 
with a spherical Gaussian density distribution. Rela- 
tively high initial electron temperature and low density 
are also required to avoid coUisional and radiative pro- 
cesses 



28, 31, 38, 3^ that would heat electrons or lead 
to electron-ion equilibration. In future studies, devia- 
tions from the model can be used to study these coUi- 
sional processes when they become important. The ex- 
pansion dynamics shown here provide a general tool for 
developing a greater understanding and intuition for the 
expansion of plasmas. 
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